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We propose a simple stochastic model of cascading transport in wave number space to clarify the 
origin of intermittent behavior of fully-developed fluid turbulence. In spite of lack of nonlinearity 
and viscosity the model gives non-Gaussian fluctuations and multifractal scalings consistent with 
experimental data. 
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Spatio-temporal structure of fully-developed isotropic 
fluid turbulence is one of frontier topics in both fluid 
mechanics and statistical physics. Although several new 
concepts, for example, fractals |2), chaos and in- 
termittency, have been introduced to characterize fully- 
developed turbulence, our present understanding of fluid 
turbulence is still far from complete. 

There are many difficulties both in experimental and 
numerical approaches. Experimentally, typical available 
data is time sequential data of velocity at a fixed point 
which is obviously not sufficient to observe global 
spatio-temporal structures of fluid turbulence. Instanta- 
neous local velocities at arbitrary points can be obtained 
by using direct numerical integration however, fully- 
developed turbulence requires so many degrees of free- 
dom that even a super computer is not powerful enough. 
To help our understanding of the behavior of turbulence 
we need simple models which are easily accessible and 
share some basic properties with fluid turbulence. 

In the history of study on turbulence one of the most 
successful approaches is Kolmogorov's scaling argument 
[01, in which the energy spectrum, E{k), was predicted 
to have a power law wave number dependence as k^^/^ . 
His prediction was based on a very simple assumption of 
cascading process that conserved energy is transported 
continuously in the wave number space toward higher 
wave numbers. Although his simple assumption captured 
the most basic property of fully-developed turbulence, a 
modification was proposed by himself and others to in- 
clude the effect of spontaneous spatial inhomogeneity so 
called the intermittency 1^ . In order to take into account 
the fluctuations of energy dissipation rates a log-normal 
distribution was assumed and the exponent of the energy 
spectrum was modified ||^. 

Another approach to the intermittency is based on ge- 
ometrical models that an eddy breaks into smaller eddies 
forming a fractal configuration [B| . This fractal approach 



chao-dyn/9503004 ,TITCMT-95-7 



is recently extended to apply the concept of multifractals 
Q and now the agreement with experiments becomes 
quantitatively acceptable jl^ 

Recently, much experimental interest is focused on 
observation of probability density function ( in short 
PDF ) of various variables in turbulence. It is well- 
confirmed that the PDF of velocity differences between 
a pair of points separated by a distance clearly deviates 
from Gaussian and the deviation becomes more evident 
for closer pairs. In terms of Fourier space this means 
that higher wave number components of velocity field 
show larger deviation from Gaussian [Q . 

The multifractal concept and the non-Gaussian PDFs 
are consistent representations of the intermittency. Actu- 
ally, by combining these two ingredients She and Leveque 
pd| have derived an analytical formula of multifractal ex- 
ponents which fits with experimental values very nicely. 
An interesting point in their derivation is that Navier- 
Stokes equation is not used explicitly like other theories 
on intermittency. 

Generally speaking lack of dynamic equation in the- 
ory apparently shows incompleteness of the theory itself, 
however, in the present case there is a positive aspect that 
the theory's applicability is not limited to Navier-Stokes 
turbulence. Actually, almost identical non-Gaussian 
PDFs and the multifractal exponents have been reported 
in an experiment on thermal convective turbulence which 
is governed by Boussinesq equation . 

Similar intermittent behavior can also be found in a 
drastically simplified numerical model called GOY model 
|]l3| , ^ . It belongs to so-called shell models as the model 
discards infinite degrees of freedom in Navier-Stokes 
equation and uses only finite numbers of representative 
wave number components called shells. It is described 
by the following set of nonlinear differential equations for 
complex velocity Fourier components, u{j, t), for discrete 
wave numbers kj — V , where 6 is a positive constant and 
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j is an integer. 

+ i {akju*{j + l,t)u*{j + 2, t) 
+bk,^iu*ij~l,t)u*{j + l,t) 

Here, the first term in right hand side shows a viscous 
dissipation effect, and the second term corresponds to 
external forcing. The remaining three terms are nonlin- 
ear terms which are designed to have similar nonlinearity 
with Navier-Stokes equation. By these nonlinear terms 
this deterministic system shows chaotic behavior and the 
resulting statistics agrees with real fluid intermittency 

In this letter, we are going to clarify why such a simpli- 
fied model can produce intermittent behavior quite sim- 
ilar to those in real turbulence having infinitely many 
degrees of freedom. Our answer is simple: The inter- 
mittency can be a very general and universal behavior 
in a wide class of nonequilibrium systems not specified 
to fluid turbulence. The key point is a fluctuating lo- 
cal directional transfer in the wave number space, so the 
number of degrees of freedom, the details of nonlinearity 
and the types of dissipation are not crucial. 

In order to validate this scenario we introduce a new 
model, the stochastic shell model, which is even much 
simpler than GOY model. Instead of the chaotic non- 
linear interactions we introduce discrete stochastic pro- 
cesses using real variable u as follows; 

■u{j,t+ 1) = u{j,t) 

+e[R{j,j-l;t)u{j-l,t)- I R{j + l,j;t) \ u{j,t)], (1) 

Here, is a non-negative constant less than unity, and 
R{j,j — is a random number which expresses the 
momentum transfer from to kj; 
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+1, Prob p/2 
-I, Prob p/2 
0, Prob 1 — p 



It should be noted that this model can be viewed as a 
multiplicative process of random matrices, thus it does 
not have any nonlinearity. Similar linear multiplica- 
tive stochastic process has been analyzed theoretically 
by Deutsch jl^ in which an anomalous scaling relation 
like the multifractal scaling is reported. 

In the following, we adopt the boundary condition that 
u{N,t) = for t > where N is the total number of 
shells. Compared with the ordinary exponential decay 
by viscous dissipation in large wave number regions, this 
artificial boundary condition of a sudden cut off may look 
too rude. However, the boundary condition at the largest 
wave number does not affect the intermediate wave num- 
ber components because our model introduces the one- 
way transport towards the larger wave numbers. Thus 



u{N,t) = is simply for numerical convenience. As 
an energy injection effect by an external forcing we add 
unifrom random numbers (e [—0.5, 0.5]) to the 0-th shell 
at every time step. 

Analytical treatment of this stochastic model can be 
done by introducing the characteristic function Yj (p, t) = 
{exp[ipu{j,t)]) where {O) denotes the averaged value of 
the variable O over realizations. Using this definition we 
get the following equation for the characteristic function 
from Eq.d). 

Y,{p,t + 1) = {(1- p)Y, {p, t) +pY,[(l- e)p, t] } 
X {{l-p)+pY,^^{ep,t)}. 

We can show analytically that the system quickly con- 
verges to a statistically steady state where all averaged 
quantities become independent of time steps. This im- 
plies that the ensemble average and a time average gives 
the same result independent of initial conditions. By ex- 
panding the characteristic function in terms of p up to 
the second order we have a rigorous steady state relation 
for 'energy', E, = (| u{j,t) ^ /2), as E, = 9E,^^/{2-e). 
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FIG. 1. Symbols: Higher moment of u{j,t).{a) q = 2 (b) 
= 6 (c) g = 10 ( x:b=16,+:b=64). Tliick lines: C<j = g/3, 
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The energy spectrum follows the familiar power law 
form, Ej = k~'^ , by choosing 9 to be 2/(1 + &*»). We em- 
ploy C = 2/3 so as to satisfy the Kolmogorov spectrum. 
The parameter p does not appear in the energy relation 
and it turns out that p is irrelevant in the following dis- 
cussion as far as < p < 1, so we fix its value to be 0.5 
hereafter. Thus, the model has only one free parameter 
b. 

To measure the intermittency quantitatively we calcu- 
late higher order moments of t) and generalize the 
definition of the scaling exponent C as (| u{j,t) I"?) oc 
fcj Figure ^ shows the results for q — 2, 6, and 10 
{N = 10, 6 = 16, 64). Here averages are taken over 2 x lO'^ 
steps. ^2 is equal to 2/3 as expected by analytical treat- 
ment. 
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FIG. 2. Comparison of PDF between stochastic shell model 
and experiment {<>). All data are normalized so as to have 
variance of unity for comparison. Top: Experiment: 0.15kHz, 
model: j — 0. Middle: Experiment: 3kHz, model: j — 1. 
Bottom: Experiment: lOkHz, model: j — 2. 

On the other hand, values of for g = 6 and 10 clearly 
deviate from Kolmogorov's original scaling exponents, 
g/3. We can conclude that our model exhibits an anoma- 
lous scaling 1 1^ . In order to compare the exponents with 
experimental values, we also plot the empirical multifrac- 
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to be a good approximation to the experimental values 
[pd| . Our data agree with this empirical values in the 
whole range, meaning that our model can reproduce the 
multifractal scaling quantitatively. 

Next, we check our model whether the wave number 
dependence of PDFs observed in experiments can also be 
reproduced correctly or not. Figure ^ shows a compar- 
ison between experimentally observed PDFs and those 
of our model (6 = 64, = 5). The experimental tur- 
bulence data are PDFs of frequency-band-pass-filtered 
velocity signals obtained by one of the authors TK [Q. 
In the case of our numerical model we plotted t) 
because the imaginary part follows the same statistics. 
Data are sampled every two steps over total 2 x 10^ steps. 
The experimental PDFs show a tendency that the devia- 
tion from Gaussian becomes more evident for larger wave 
numbers. The numerical data fit them very nicely in the 
whole range in each case. Namely, our model also repro- 
duces the other aspect of intermittency, the non-Gaussian 
PDFs @. 

In high wave number limit it is known that both 
GOY model and Navier-Stokes turbulence have almost 
identical power laws in PDF P{u) of velocity u, i.e., 
P(u) ~ with (3 ~ 1.6 ^7^. Thus, the appearance 
of power tails in PDFs at very high wavenumber may 
also be an important nature of fluid intermittency. 




tal relation, (q — q/9 + 2 



1 



(I) 



which is known 



FIG. 3. Power law PDF obtained by the stochastic shell 
model {j = 40). The straight line indicates P{u) ~ 

Figure ^ shows PDF of our model in very high wave 
number region {j — 40, N — 50, b — 16). Data are sam- 
pled every two steps over total 2 x 10^ steps. The PDF 
shows a power law with its exponent about —2. Taking 
into account the errors accompanied with the estimation 
of such delicate exponent, this value is not far from the 
values of fluid turbulence, —1.6. We judge that our model 
is also consistent on this power law behavior jl^]. Re- 
membering that our model is linear and has no viscosity, 
the above results suggests that the intermittency may 
not have any direct relationship with nonlinearity and 
viscosity. This may sounds quite contrary to the com- 
mon sense, and a big question naturally arises, that is, 
what is really responsible for the intermittency? 

To elucidate minimal factors for the realization of in- 
termittency we still modify the model. Instead of real 
variables we consider complex variables for u{j, t) in our 
stochastic shell model. The results are almost identical 
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with what we have reported. However, when we fix the 
sign of R{j,j — 1; f); independent of whether u is real or 
complex, we get considerably different state. We have a 
statistically steady result and the distributions deviate 
more from Gaussian for larger wave number components 
as in the case with sign randomization, but the energy 
spectrum is modified to a non-power law and the mul- 
tifractal scalings are lost completely. This difference is 
expected to be caused by the strong correlations among 
velocity components in each shell. This result demon- 
strates that transport with sign (or phase for complex u) 
randomization is essential for the intermittency in real 
turbulence. 

Concluding the paper, we have introduced a stochas- 
tic shell modcJ to describe cascading transports in wave 
number space, which can be viewed as a kind of lin- 
ear multiplicative stochastic process. The model with 
sign randomization reproduces the intermittent behavior 
correctly while its constant sign version looses all such 
behavior. This fact implies that the origin of intermit- 
tency is in the cascading process with sign (or phase for 
complex u) . Nonlinearity and viscosity does not ap- 
pear explicitly in our discussion, therefore, these effects 
may give only indirect contribution to the intermittency. 
Non-linearity in real turbulence is important as a origin of 
randomness. However, once randomness is introduced, it 
is not important whether randomness originates in non- 
linearity or white noise when we consider stochastic prop- 
erty like intermittency and anomalous scaling. Further 
investigation of the stochastic shell model will contribute 
on the problem of universality of intermittency. 
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